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Abstract. We study the continuity of an abstract generalization 
of the maximum-entropy inference — a maximizer. It is defined as 
a right-inverse of a linear map restricted to a convex body which 
uniquely maximizes on each fiber of the linear map a continuous 
function on the convex body. Using convex geometry we prove, 
amongst others, the existence of discontinuities of the maximizer 
at limits of extremal points not being extremal points themselves 
and apply the result to quantum correlations. Further, we use 
numerical range methods in the case of quantum inference which 
refers to two observables. One result is a complete characterization 
of points of discontinuity for 3 x 3 matrices. 
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1. Introduction 

The maximum-entropy principle is a topic in physics since the 19 th 
century in the work of Boltzmann, Gibbs, von Neumann, Jaynes and 
many others [40, 23]. The continuity issue of the maximum-entropy 
inference under linear constraints in the quantum context was studied 
by one of the authors in [48] (jointly with A. Knauf), and then further 
in [44]—[46]. Discontinuity points of the maximum-entropy inference 
are a distinguished quantum phenomenon because they are missing in 
the analogous maximum-entropy inference of probability distributions 
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which is formally included in the quantum setting by viewing probabil¬ 
ity vectors as diagonal matrices. The discontinuity points were recently 
discussed by Chen et al. [13] in condensed matter physics. 

Convex geometry has proved powerful for the problem of disconti¬ 
nuity in our preceding contribution [45]. The crucial difference with 
the simplex of probability distributions is that a quantum state space, 
consisting of density matrices, has curved and flat boundary portions. 
Already its planar linear images have a very rich geometry. They corre¬ 
spond to the notion of numerical range in operator theory whose shapes 
are well understood for 2 x 2 (the classical elliptical range theorem) and 
3x3 matrices, see the work of Kippenhahn [24] and earlier papers by 
the first two authors [26] (jointly with D. Keeler) and [36]. 

Let us briefly return to the condensed matter note which is a sub¬ 
stantial motivation for this work. Definitions of a phase transition in 
statistical mechanics include symmetry breaking, such as melting crys¬ 
tals, and long-range correlations which may be certified by a power law 
correlation function, see for example Yeomans [51]. Quantum phase 
transitions are not necessarily associated with symmetry breaking or 
long-range correlations, for example see the discussion by Wen [49], 
Sec. 1.4. Still, correlations explain quantum phase transitions in some 
cases: Liu et al. [31] have recognized that a quantum phase transition 
in Kitaev’s toric code model can be seen in the six-body correlations of 
the ground states of a 3 x 4 torus model. We demonstrate by an exam¬ 
ple that our methods are suitable to obtain analytical results about the 
correlation quantities used by Liu et al. which are known as irreducible 
many-party correlations and which were defined by Linden et al. and 
Zhou [30, 52] based on the maximum-entropy principle. 

The aim of this article is to contribute to the continuity theory of 
the maximum-entropy inference and of the irreducible correlation us¬ 
ing techniques from convex geometry and operator theory (numerical 
range). Convex geometrical methods will be developed in reference to 
a maximizer H which is defined as a right-inverse of a linear map / 
restricted to a convex body K which uniquely maximizes on each fiber 
of ./1 K a continuous function on K. We have in mind the example 
where H is an abstract generalization of the maximum-entropy infer¬ 
ence (defined in Sec. 2) and the continuous function is a generalization 
of the von Neumann entropy. 

Our analysis will be based on our continuity result [45] which we 
recall in Sec. 2 and which allows us to study the continuity of the 
maximizer without solving the respective inverse problem explicitly. 
This is possible by studying the openness of the restricted linear map 
f\i<- Continuity results follow as corollaries of openness results. 
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The results presented in Sec. 3 have appeared earlier in the paper 
[45] by the fourth author and are based on gauge functions of the 
domain L := f(K) of the maximizer. Here we provide a new unified 
proof in terms of the notion of simplicial point which is a point-wise 
defined variant of a locally simplicial set in the sense of Rockafellar 
[35]. We prove that the maximizer is continuous at all simplicial points 
- in particular it is continuous in the restriction to a polytope or to a 
relatively open convex subset. 

Further, we present two new results inspired by examples by Chen 
et al. [13]. In Sec. 4 we prove a dichotomy with regard to the partition 
of the domain L of the maximizer II : the continuity of H is equivalent 
to the continuity of its restriction to the relative boundary of L and 
can be decided in terms of its restriction to the relative interior of L. 

In Sec. 5 we prove a necessary condition for the continuity of H 
in terms of the face function studied by Klee and Martin [25] and 
others. The face function maps every point of L to the unique face of 
L which contains the given point in its relative interior. For example, 
we prove under mild assumptions (which are satisfied by the maximum- 
entropy inference) that the maximizer is discontinuous at points w G L 
which are limit points of extremal points of L but not extremal points 
themselves. In that case we remark in Exa. 5.2 that the discontinuity at 
w is not removable from the restriction of H to the relative boundary of 
L by changing only the value at w. Discontinuities of the (unrestricted) 
maximizer H are never removable as we will point out in Sec. 4. 

In the three last sections we specialize to the convex body K equal 
to a quantum state space. The aim of Sec. 6 is to demonstrate that 
convex geometry is an essential part of the topology of quantum corre¬ 
lations. We use the face function method to point out discontinuities 
in the irreducible three-party correlation of three qnbits which is possi¬ 
ble because this correlation quantity is based on the maximum-entropy 
principle. We check onr analysis by consulting a result in the context 
of pure state reconstruction by Linden et al. [30]. 

In Secs. 7 and 8 we consider two hermitian d x d matrices, del, 
having the meaning of quantum observables. Then the domain L of the 
maximizer is planar and can be identified with the numerical range of a 
complex d.xd matrix A, that is W(A) := {x*Ax \ x G C S d } where C S d 
denotes the unit-sphere in C d . We determine the maximum number of 
discontinuities in dimension d — 4,5 (new methods are needed for 
d > 6). Further, in Sec. 8 we give a complete characterization of points 
of discontinuity of the maximum-entropy inference for d — 3. These 
results essentially are corollaries of the Theorem 7.3 which connects to 
the paper [28] by T. Leake, B. Lins and the second author, devoted 
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to the numerical range. One of the proof ideas is that the notion of 
simplicial point, mentioned above, is a complement of the notion of 
round boundary point in the second author’s papers [15] (jointly with 
D. Corey, C.R. Johnson, R. Kirk and B. Lins) and [27, 28] (jointly 
with T. Leake and B. Lins) on the numerical range. The second idea is 
concerned with singularly generated points which were defined in [28] 
as having exactly one linearly independent unit vector in the pre-image 
under the map C S d —> C, x (->• x*Ax. 

It is worth mentioning that similar pre-image problems of the nu¬ 
merical range are of a broader interest, see for example Carden [9]. In 
particular, pre-image problems for more than two observables appear in 
quantum state reconstruction, see for example Gross et al., Heinosaari 
et al. and Chen et al. [17, 19, 12], and in quantum chemistry, see for 
example Erdahl, Ocko et al. and Chen et al. [16, 33, 11]. 

Acknowledgements. AS and SW appreciate the support by the German Research 
Foundation within the project “Quantum statistics: decision problems and entropic 
functionals on state spaces” (10/T1-09/T4). IS was supported in part by the 
Plumeri Award for Faculty Excellence from the College of William and Mary and 
by Faculty Research funding from the Division of Science and Mathematics, New 
York University Abu Dhabi. 


2. Preliminaries 

In this section we introduce the central objects and concepts of this 
article including the notions of state space and maximum-entropy infer¬ 
ence in quantum mechanics, starting with their abstract counterparts 
of convex body and maximizer, as well as the joint numerical range 
and the topological notion of open map. 

Definition 2.1 (Maximizer). Let X, Y be finite-dimensional real norm- 
ed vector spaces, let / : X —> Y be a continuous map and let K c X 
be compact. The compact set L := f(K) C Y parametrizes the fibers 

f\lc( w ) = i v e K I f( v ) = w e 

of f\i<- We assume the inverse problem of selecting a point in each 
fiber is solved by maximizing a continuous function g : K — > R which 
attains a unique maximum on each fiber. The maximizer is defined by 

H : L —> K, w i— y argmax{g(u) \ v E f\j/ H}. 

All sets which will be introduced in the sequel are tacitly assumed to 
be subsets of a finite-dimensional real normed vector space. Unless 
otherwise stated we will always assume that / is linear and that K is 
a convex body , that is a compact and convex set. 
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We remark that the domain of / is X rather than K by consistency 
with our main example (2.1). Since K is compact, this choice is no 
restriction because, by the Tietze extension theorem, any continuous 
function K —y Y can be extended to a continuous function X —y Y. 

On K c X and L = f(K) C Y we use the subspace topology 
induced by X and of Y, respectively. A neighborhood of a point x 
in a topological space X is any subset of X containing an open set 
containing x. We call a map 7 : X —>■ Y between topological spaces 
X, Y open at x G X if for any neighborhood N C X of x the image 
7 (N) is a neighborhood of -/(x) in Y. We say 7 is open on a given 
subset of X if 7 is open at each point in the subset. Finally, 7 is open 
if 7 is open on X. 

In Thm. 4.9 in [45] one of the authors has proved the following. 

Fact 2.2 (Continuity-openness equivalence). Let K C X be an arbi¬ 
trary compact subset, not necessarily convex, and let / : X —>• Y be 
an arbitrary continuous function, not necessarily linear. Then for any 
w G L the maximizer H is continuous at w if and only if f\x is open 
at H(w ). 

From now on we make the global assumption that K is a convex 
body and that / is linear. We will argue in terms of the openness 
of f\i<- All results may, and some will, be translated into continuity 
statements of H using Fact 2.2. 

Our main example of convex body K will be the state space of a 
matrix algebra which is studied in operator theory [1], Let M dj d G N, 
denote the full matrix algebra of d x d-matrices with complex coeffi¬ 
cients. The algebra M d is a complex C*-algebra with identity 1^. We 
shall also write 0 = CL for the zero in M d . Let A denote a (complex) 
C*-subalgebra of M d . For example, we will introduce in Fact 5.1 the 
C*-algebras pM d p = {pap \ a G M d } where p G M d is a projection, 
that is a hermitian idempotent p = p* = p 2 . See Lemma 9.1 for other 
relevant examples of C*-algebras. 

We denote by A h — {a G A \ a* — a} the real vector space of 
hermitian matrices in A, known as observables in physics, and we endow 
it with the scalar product (a, b) := tr (ab), a,b G A h which makes A h a 
Euclidean space. We call state space of A the convex body 

M(A) := {p G A | p Y 0, tr(p) = 1}. 

Here a >z 0 means the matrix a G A is positive semi-definite, that is 
a* = a and all eigenvalues of a are non-negative. Elements of M.(A) 
are called density matrices in physics [1, 5, 32], They are in one-to-one 
correspondence to the positive normalized linear functionals A —» C 
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called states, see for example [1], Sec. 4. We use the terms state and 
density matrix synonymously. 

Now we confine definitions to the full matrix algebra M d and we write 
Aid '■= M{M d ). Given a number r G N of fixed observables Ui G M§, 
i = 1 ,..., r, we write u = [u\,..., u r ) and we define the expected value 
function 

(2.1) E = E u : Mj — y M r , a (->• ((rq, a), ..., (u r , a)). 

In our earlier papers [42, 45] we have called the set of expected values 

L( u) = L(u u ...,u r ) := (E u (p) | p G M d } C M r 

the convex support. This name is motivated by probability theory [3]. 
The probability vectors of length d, embedded as diagonal matrices 
into M d , are the states of the algebra of diagonal matrices. Random 
variables on {1 ,... ,d} correspond to diagonal matrices U\,... ,u r and 
the set of their expected value tuples, called convex support in [3], 
equals L(u) by Lemma 9.1. 

Given expected values a G L(u) the maximum-entropy state p*(a) 
is the unique state in the fiber E|j^ (a) which maximizes on E|^ (a) 
the von Neumann entropy 

S(p) = -tr(p-lnp). 

Functional calculus with respect to the continuous function [0,1] —y M, 
x i— y x ' ln(x) where 0 ■ InO = 0 can be used to define S. The mapping 

(2.2) p* L(u) —y Midi Q! i —y p*{of) 

is called the maximum-entropy inference, see [23, 50, 22] for more de¬ 
tails. In physics, the von Neumann entropy quantifies the uncertainty 
in a state [41], The state p*(a) is considered the most non-committal, 
most unbiased or least informative state with regard to all missing 
information beyond the expected values a [23]. 

Although p* can be discontinuous it is smooth up to boundary points. 
We denote by C the norm closure of any set C. The relative interior 
ri(C') of C is the interior of C in the topology of the affine hull of C, 
and C is relatively open if C = ri(C) holds. The relative boundary of 
C is rb(C) := C \ ri(C). The following statement is proved in [50], 
Thm. 2b. 

Fact 2.3 (Real analyticity). The maximum-entropy inference p*(a) is 
real analytic in the relative interior of L(u). 

Now we turn to the joint numerical range which will be useful to 
address the continuity of p* for two (r = 2) observables. It will be 
convenient to denote the inner product of x,y G C d by x*y := aqq/i + 
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• • • + XdVd and to denote by xy* : C d —> C d the linear map defined for 
z E C d by (xy*)(z) := ( y*z)x . Vectors in C n , n E N, will be understood 
as column vectors. To save space we will write them equivalently in 
the column and row forms 

(xi,...,x n ), Xi,... ,x n E C. 


Definition 2.4. The joint numerical range of u = (rq,... ,u r ) is the 
subset of M r defined by 

W (u) = W(u \,..., u r ) : = {( x*U\X ,..., x*u r x) \ x E C d , x*x = 1}. 

Let us recall that the convex hull of the joint numerical range is the 
convex support. We denote the convex hull of any set C by conv(C). 
For all d,, r eN we have 

(2.3) conv(VF(u)) = L(u). 

See [18], Thm. 1, for the identity (2.3) formulated as an affine isomor¬ 
phism. For two observables (r = 2), if we identify M 2 = C, then the 
joint numerical range W (iq, m 2 ) equals the numerical range W (tq + i-iq) 
which is convex by the Toeplitz-Hausdorff theorem. Hence (2.3) implies 

(2.4) L(ui,u 2 ) = W(ui, u 2 ). 

A proof of (2.4) can be found in [6], Thm. 3. For three observables 
(r = 3) the joint numerical range W(ui,u 2 ,u 3 ) is also convex but only 
for matrix size d > 3, see [2, 29]. 

Let us discuss easy properties and example of the convex support. 

Remark 2.5. The following transformations do not essentially alter 
L(u), W(u), and p*{ot). 

(1) Remove any u t which is a (real) linear combination of iq, ..., 

U{ + 1 , ..., u r . Conversely, add observables which are linear combi¬ 
nations of tq,..., u r . 

(2) Replace Ui with + q 1, where c* E M, for any i. 

(3) Replace Ui with T*u{T for all i, where T E M c i is a unitary. 

The state space JA 2 of M 2 is a three-dimensional Euclidean ball, 
known as Bloch ball [5]. Its surface is known as the Bloch sphere [5, 32] 
or Poincare sphere [1]. The openness of E|cpi on the Bloch sphere is 
proved in Coro. 6 in [15]. The openness of E|_a/< 2 on the Bloch ball is 
shown in Example 4.15.2 in [45]: 


Xi 

x 2 
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Fact 2.6. The expected value function E |^ 2 is open. 


The possible convex support sets L(u) of Xi 2 are the linear im¬ 
ages of the Bloch ball and they are easily identified algebraically. By 
Rem. 2.5(1,2) we may assume that r < 3 and that ui,u 2 ,u 3 G M £ are 
zero-trace and mutually orthogonal. If r = 3, then by Rem. 2.5(3) we 
can take the observables equal to the Pauli matrices 



' 0 

1 ' 


' 0 ' 


n 0 1 

: = 

1 

0 

, ^2 : = 

i 0 

, 03 : = 

0 

1 


It is well known (and easy to check) that L(cr 1 , a 2 , cr 3 ) is the unit ball 
in M 3 . Notice that the joint numerical range W(ai, cr 2 , 03) is the unit 
sphere in M 3 which is not convex. If r = 2, then we can take U\ = 07 , 
u 2 = <j 2 . Then L( u) is the numerical range of 


07 + ia 2 


0 2 
0 0 


which is the unit disk in C centered at zero. If r = 1, then L(u) is a 
line segment (we ignore the trivial case when L(u) is a singleton). 


3. Simplicial Points 

In this section we recapitulate results from [45] now with a unified 
proof in terms of simplicial points. The term of simplicial point has 
another advantage that it complements the term of round boundary 
point in Lemma 7.1. 

The main idea is the gauge condition in Fact 3.1 from our work [45] 
which is somewhat similar to the idea of our Thm. 4 in [15]. Let X be 
a finite-dimensional real norrned vector space and let C be a convex 
subset of X. The gauge of C is defined by 

7 c(v) := inf{A > 0 | v G AC}, v G X. 

Recall that the gauge of the unit ball in X is the norm. More generally, 
7 c is positively homogeneous of degree one and convex [35]. If C 7 ^ 0 
then the positive hull of C is defined by pos(C) := (An | A > 0, v G C}. 

In [45], Prop. 4.11, one of the authors has proved the following. 

Fact 3.1 (Gauge condition). If w G L and if the gauge 7 l-w is bounded 
on the set of unit vectors in pos(L — w) with respect to an arbitrary 
norm on Y then /1 k is open on f\j/(w). 

The assumptions of Fact 3.1 are fulfilled for example at the apices of 
the skew cone (5.1). The following Prop. 3.2 does not apply there but 
it will suffice for the purposes of this article. 
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We call a point x in a convex set C a simplicial point if there exists a 
finite set of simplices Si,, S m C C such that the union Si U... U S m 
is a neighborhood of x in C. A convex set C is locally simplicial [35] if 
all its elements are simplicial points. 

Proposition 3.2 (Simplicial points). If w G L is a simplicial point of 
L then f \k is open on f\f c 1 (w). 

Proof: Let Si,..., S m C L, m E N, be a set of simplices such that 
the union U Si U ... U S m is a neighborhood of w in L. Since L 
is convex, the convex hull P of U is also a neighborhood of w in L. 
Therefore pos(L — w) = pos(P — w ) holds. Since P — w C L — w holds 
we have for all vectors u G pos(P — w ) the inequality 

7 L-wiu) < 7 p- w {u). 

For unit vectors u in pos (P—w) the right-hand side is bounded because 
P — w is polyhedral convex and contains the origin, see Rem. 3.1 in 
[39]. Therefore Fact 3.1 implies the claim. □ 

We mention some examples where we will apply Prop. 3.2. The 
relative interior of L and polytopes included in L are locally simplicial 
sets [35]. So the maximizer H is continuous on the relative interior of 
L and globally continuous if L is a polytope. Moreover, the restriction 
H\p is continuous for every polytope P C L. For example, the convex 
support L(u) is a polytope for commutative observables ui,...,u r , see 
Sec. 2 in [45]. 

As the last example we mention that relative interior points of facets 
are simplicial points and we leave the proof to the reader because the 
openness of f\x on the fibers of these points is also proved in Coro. 4.4 
in the next section. 

Recall that a face [35] of a convex set C is any convex subset F C 
C which contains all segments in C which meet F with an interior 
point. A face which is a singleton is called an extremal point. A face 
of codimension one in C is a facet of C. 

4. Boundary-Interior Dichotomy 

We show that the continuity of the maximizer H is certified by its 
restriction to the relative interior ri(L) of L and also by the restriction 
to the relative boundary rb(L) of L. 

We begin with the relative interior by citing from Lemma 4.8 in our 
paper [45]: 

Fact 4.1 (Norm closure). We have H(L) C H(n(L)). 
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This statement was proved earlier in the context of the maximum- 
entropy inference (2.2) in Thm. 2d in [50]. The continuity of the max¬ 
imizer H can be decided using H | ri ( L ). In fact, given w G L, if 

is a singleton {v} then H is continuous at w and H{w ) = v. Otherwise 
H is discontinuous at w. Fact 4.1 shows also that discontinuities of H 
are not removable. 

Turning to the relative boundary we will use the property that for 
every w G rb(L) the convex hull of a neighborhood of w in rb(L) and 
of a point in ri(L) is a neighborhood of w in L. This is easy to check 
for a Euclidean ball L with center w. The following fact, proved in 
Sec. 8.1 in [7], generalizes this from the ball to arbitrary convex bodies. 
Recall from [35] that a mapping 7 : M n —> W 1 , n G N, is positively 
homogeneous of degree one if for each x G M n we have 'y(Xx) = Ay/c), 
0 < A < 00 . 

Fact 4.2 (Thm. of Sz. Nagy). Let C C M n , n G N, be a convex body 
containing the origin in its interior. Then there exists a homeomor- 
phism from C onto the standard Euclidean unit ball of which is 
positively homogeneous of degree one. 

Theorem 4.3. Let K := /1 ff ( rb(L)). For all v G K the map f \k is 
open at v if and only if f \^ is open at v. 

Proof: One direction follows by taking intersections of neighbor¬ 
hoods. Let us prove conversely that the openness of f\^ at v G K 
implies the openness of / \k at v. 

Let N C K be a neighborhood of v in K. Then N 0 K is a neigh¬ 
borhood of v in K. By assumptions f\^ is open at v so the image 
f(N fi K) is a neighborhood of f(v) in rb(L). Since X is locally con¬ 
vex, we can assume that N is convex so conv(v, N D K) C A r for a 
given point v G N \ K. The linearity of / shows 

(4.1) /(IV) D /(conv(u, N D K)) = conv(/(u), f(N fl K)). 

Since f(v) lies in the relative interior of L the discussion in the para¬ 
graph before Fact 4.2 proves that conv(/(v), f(N fl K)) is a neighbor¬ 
hood of f(v) in L. Then (4.1) shows that f(N) is a neighborhood of 
f(y) in L which completes the proof. □ 

Thm. 4.3 applies to any facet F of L because F is a neighborhood 
in rb(L) of the relative interior points of F. 
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Corollary 4.4. Let F be a facet of L and let w G ri (F). Then f \ k is 
open on f\j/(w). 


5. The Face Function of L 

We prove a necessary continuity condition of the maximizer H in 
terms of the lower semi-continuity of the face function of L. The lower 
semi-continuity of the face function implies that a limit of extremal 
points is an extremal point [34], 

An example where a limit of extremal points is not an extremal point 
is given by the convex hull of 

(5.1) {(s, t, 0) G M 3 | (s - l ) 2 + t 2 = 1} U {(0, 0, ±1)}. 

Here the set of extremal points {(s, t, 0) Gl 3 | (s - l ) 2 + t 2 = l,s^0} 
contains ( 0 , 0 , 0 ) which is an interior point of the segment connecting 
( 0 , 0 ,- 1 ) and ( 0 , 0 , 1 ). 

To provide an example with the state space M.% we need algebraic 
representations of faces. A subset F of a convex set C is an exposed 
face of C if F = 0 or if F equals the set of maximizers in C of a linear 
functional. One can show that every exposed face of C is a face of C. 
An exposed extremal point is called exposed point. 

Fact 5.1 (Faces of state spaces). 

(1) The non-empty faces of the state space Add are of the form A4(pMdp) 
where p G Aid is a non-zero projection, see for example [42], Sec. 2.3. 

(2) Consider the exposed face F = argmax{a;*A | a G L(u)} of L(u) 
where A = (Ai,..., A r ) G M r . Then E|j^ d (F) = A4(pM d p) where p 
is the spectral projection of the matrix u(A) = Ai u\ + • • • + A r u r 
corresponding to the maximal eigenvalue of u(A). This follows from 
(p, u(A)) = E(p)*A, p G A4 d , A G M r , and from [42], Thm. 2.9. 

(3) Extremal points of the state space Ai d are called pure states. A 
state is a pure state if and only if it is of the form xx* for some unit 
vector x G C d , see for example [1], Sec. 4. 

We now discuss a three-dimensional linear image of Ai% which has 
appeared as Exa. 4 in [13]. This linear image has a sequence of extremal 
points which converge to a point which is not an extremal point. A 
three-dimensional cross-section of .M 3 with this property is discussed 
in Rem 5.9 in [45]. 
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Example 5.2. We consider the convex support L(ui, U 2 ,u 3 ) of 


Ui := 

'10 0 

0 1 0 

, U 2 ■■ = 

'10 1 

0 1 1 

, u 3 := 

'10 1 

0 0 1 


0 0-1 


11-1 


11-1 


A limit of extremal points. Let e G M and £(e) := y/Jl — e) 2 + 2e 2 . 
The eigenvalues of r u.\ — eu 2 are {1 — e, ±£(e)}. If e 0 then the 
maximal eigenvalue £(e) is non-degenerate and v(e) := (1,1, e • x(e)) 
is a corresponding eigenvector, where x(e) := (£(e) + e — l)/e 2 . By 
Fact 5.1(2) the pure state p(e) := c(e) • v(e)v(e)*, where c(e) > 0 is for 
normalization, defines the exposed point 

a(e) := E (p(e)) = 2 _ff )3!(e) (1 - e, 1 - 3e, 1 - 3e) - (0,0, 2+(e .* (e))a ) 
of L(u). Since x(e) —>■ 1 for e —>■ 0 we have 

ct( 0 ) := liirif —> 0 ot(e) = ( 1 , 1 , |). 

We could also do the easier computation a(0) = E(lim e _ > . 0 p(e)) but the 
focus should be on L(u) rather than A4 3 . The limit a(0) is not an 
extremal point because it is the mid-point of the segment s C L(u) 
between 

E((0,1,0)(0,1, 0)*) = ( 1 , 1 , 0 ) and E((l, 0,0)(1, 0, 0)*) = (1,1,1). 

Discontinuous maximum-entropy inference. We would like to point 
out that the maximum-entropy inference p* is discontinuous at a( 0 ) 
along the curve a(e). Namely, for e ^ 0 we have p*(a(e)) = p(e) by 
Fact 5.1(2). The limit 

linv^o p*(a(e)) = |(1,1,0)(1,1, 0)* 

is a pure state while p*(a(0)) has rank two. Indeed, we obtain 

p>( 0 )) = p /2 for p := ( 0 , 1 , 0 )( 0 , 1 , 0 )* + ( 1 , 0 , 0 )( 1 , 0 , 0 )L 

To see this, observe that the spectral projection corresponding to the 
maximal eigenvalue of U\ equals p. Hence, by Fact 5.1(2), the ex¬ 
posed face of L(u) consisting of the maximizers of the linear functional 
L( u) -)• I, A G (1,0, 0)*A has the pre-image M.(pM :i p) under E|_a 4 3 . 
Since a(0) = E(p/2) and since S (p/2) = log(2) is the maximal value of 
the von Neumann entropy on AA^pM^p) we have p*(a(0)) = p/2. 

Non-removability of the discontinuity. Further, we would like to 
point out that, by Prop. 3.2, the restricted maximum-entropy inference 
p*\ s is continuous on the segment s at a(0). This, together with the 
discontinuity at a( 0 ) along the curve a proves that the discontinuity at 
a( 0 ) is not removable from the restriction of p* to the relative boundary 
of L(u) (by changing only the value at a( 0 )). 
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We will show that the discontinuity in Exa. 5.2 is a consequence of 
the fact that the extremal points a(e) converge to the point a(0) which 
lies in the relative interior of a higher-dimensional face. A convex body 
C is stable [34, 14] if the mid-point map 

(5.2) CxC-iC, (x,y) ^ \{x + y) 

is open. The state space Aid is stable. Indeed, Lemma 3 in [38] proves 
that the map Aid x A id x [0,1], (p, a, A) (->• (1 — A )p + A o is open. 
Prop. 1.1 in [14] then shows that Aid is stable. 

Let us recall an equivalent statement of stability. Fact 5.3 is proved 
for example in Thm. 18.2 in [35]. 

Fact 5.3 (Face function). For each point x in a convex set C there 
exists a unique face of C which contains x in the relative interior. 

If a; is a point in a convex set C then we denote the face of C 
containing x in its relative interior simply by F(x ), omitting C (which 
should be clear from the context). The face function of C is the set- 
valued map 

C —> C, x (->• F(x) 

which has been studied for example in [25, 34], A set-valued map T : 
X —» Y between topological spaces A", Y is lower semi-continuous at 
x G X if for each open set G meeting T(x) there exists a neighborhood 
N of x such that for all 6 IV we have G fl T(a; / ) ^ 0. The set-valued 
function T is lower semi-continuous if T is lower semi-continuous at 
every point of X. 

Fact 5.4 (Stable convex bodies). If C is a convex body then the fol¬ 
lowing are equivalent [34]: 

(1) The convex body C is stable. 

(2) The face function x f—> F(x) of C is lower semi-continuous. 

(3) The function C —» N 0 , x e-X dim(F(x)) is lower semi-continuous. 

Since extremal points have dimension zero, Fact 5.4(1) and (3) prove 
for any stable convex body that a limit of extremal points must be an 
extremal point. We now list some basic relations between the face 
functions of the convex bodies K and L = f(K). 

Lemma 5.5 (Linear images of faces). Let w G L. Then: 

(1) Ifv G f\fl{w) then f{F(v)) C F(w). 

(2) Ifve ri(/| k( w )) then f ( F ( v )) = F ( w )- 
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Proof: We prove (1) assuming v G The point v is a relative 

interior point of F(v) by the definition of the face function. The relative 
interior of the linear image of a convex set is the linear image of the 
relative interior by [35], Thm. 6.6. So the relative interior of f(F(v)) is 
f(n(F(v)). Hence w = f(v) is a relative interior point of the convex set 
f(F(v)). Therefore, and since w lies in the face F(w) of L, Thm. 18.1 
in [35] proves (1). 

We prove (2) assuming v G Yi(f\ff(w)). Recall that inverse images 
of faces are faces. So G := f\ff(F(w)) is a face of K. As we have 
recalled in the previous paragraph, /(ri(G)) = ri(/(G)) holds so the 
affine space f~ 1 (w) meets the relative interior of G. Hence, by [35], 
Coro. 6.5.1, the relative interior of f\f x (w) = / -1 ( w ) fl G is the inter¬ 
section of f~ l {w) with ri(G). This shows that v lies in ri(G) and proves 
G = F(v) which completes the proof. □ 

We are ready for the main result of this section. 

Theorem 5.6. Let (wf)^ C L converge to a point w G L. We assume 
that (a) the convex body K is stable. We also assume that (b) a sequence 
(^i)ieN C K converges to a point v G K such that E (vf) = Wi, j G N, 
and such thatv lies inri(f\jf(w)). Then the following statements hold. 

(1) For all open subsets O C L meeting F(w) there exists IVgN such 
that for all i > N we have O D F{wf) 0. 

(2) We have dimF(u;) < lim inf^oo dim F{wf). 

Proof: We prove (1) assuming O C L is an open set meeting F(w). 
Since f\x is continuous O := f\f ] (O) is open. By the assumption (b) 
the point v lies in the relative interior of f\j/(w) so f(F(v)) D F(w) 
holds by Lemma 5.5(2) and O fl F(v) 0 follows. By assumption (a) 
the convex body K is stable so the face function of K is lower semi- 
continuous by Fact 5.4(1) and (2). As v = lim^oo Vi holds there is 
]VgN such that for i > N we have ODF^vf) 0. As f(F(vi )) C F{wf) 
holds by Lemma 5.5(1) we get OGF[wi) 0 for i > N. The statement 
(2) is an easy corollary of (1) for arbitrary convex bodies L. □ 

Let us discuss Thm. 5.6. 

Remark 5.7. (1) Thm. 5.6 allows us to detect discontinuities of the 
maximum-entropy inference p* in terms of the convex geometry of 
L(u). We have seen in the paragraph of (5.2) that the state space 
Aid is stable while p*(a) G ri E|[^ (a) holds for all a G L(u) by 
Lemma 5.8 in [45]. Thus the assumptions of Thm. 5.6 are fulfilled. 
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The Example 5.2 demonstrates explicitly a discontinuity of p* at 
a point a( 0 ) which is a limit of extremal points a(e) but not an 
extremal point itself. 

(2) All assumptions of Thm. 5.6 are needed in general. Consider a 
convex body K C X which is not stable, take X = Y and the 
identity map / : X —> Y. Then the face function of L = K is not 
lower semi-continuous by Fact 5.4(1) and (2), thereby contradicting 
Thm. 5.6(1). The assumption (b) is not met in Exa. 5.2 where the 
limit point p(0) = ^(1,1,0)(1,1, 0)* of p{e) G Xi 3 for e —» 0 is 
a relative boundary point of the Bloch ball E|( v |, j (o;(0)). In this 
example we have a jump from the dimension zero of F(a(e)), e ^ 0, 
to the dimension one of E(o;(0)) which contradicts Thm. 5.6(2). 

(3) Thm. 5.6(1) is in general stronger than Thm. 5.6(2) if the assump¬ 
tions (a) and (b) of the theorem are ignored. This was pointed out 
in the work [34] in the example recalled in (5.1) where (1) fails at 
all points but the vertices of the segment between ( 0 , 0 ,— 1 ) and 
( 0 , 0 , 1 ) while ( 2 ) holds at all points but the mid-point of this seg¬ 
ment. However we do not know whether the convex body in (5.1) 
is the linear image of a stable convex body so ( 1 ) and ( 2 ) could be 
equivalent under the assumptions of the theorem. 

6. Irreducible Correlation 

In this section we show that the discontinuity of the three-party 
irreducible correlation of three qubits, the correlation that can not be 
observed in two-party subsystems [30], can be detected via the face 
function of the convex body of two-party marginals. 

The discontinuity which we ’detect’ follows also from the well-known 
result [30] that almost every pure state of three qubits is uniquely 
specified by its two-party marginals among all states (pure or mixed). 
The only exceptions are the GHZ-like states [32] 

(6.1) vj) := a|000) + /5|111), |a | 2 + |/ 3| 2 = 1 , a , f3 G C 

and their local unitary transforms, that is vectors (U\ <8> U 2 8) U 3 )^ 
where U\,U 2 ,U 3 G M 2 are unitaries. 

A three-qubit system ABC is described by the algebra M 8 = M 2 < 8 > 
M 2 8> M 2 with state space M. 8 = J\4(M 2 < 8 > M 2 < 8 > M 2 ). A two-local 
Hamiltonian is a sum of tensor product terms a 8) b 8) c with at most 
two non-scalar factors a, b, c G . In this section, we fix any spanning 
set Mi,,.,, u r G Mg, r G N, of the space of two-local Hamiltonians 

and put u = (rq,..., u r ). Given any three-qubit state p G Xi 8 , its 
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marginal pab € M.\ on the AB subsystem is defined by 

{pAB,a®b) = 1 2 ), a,b G M 2 . 

The marginals pac , Pbc are defined similarly. We denote by p A = 
(Pab, Pac, Pbc) the vector of two-party marginals of p. For two states 
p, o G A4s the expected values of two-local Hamiltonians satisfy E u (<r) = 
E u (p) if and only if the two-party marginals coincide, that is a 1 ' 2 ' 1 = /h 2 d 
Thus we identify the convex support L(u) with respect to the spanning 
set u\,... ,u r of and the set of two-party marginals, that is 

(6.2) L{ u) ^ {p (2) | p G Af 8 }. 

Using the identification (6.2) we define, as in (2.2), the maximum- 
entropy inference 

p* : L( u) —>■ A4 8 , a e-)- argmax{5'(p) G A4 8 | P ^ = «}• 

The irreducible three-party correlation [30] of p is defined as the differ¬ 
ence of von Neumann entropies 

(6.3) C 3 (p) := S(p‘(p m )) - S(p). 

The definition of C 3 (p) is derived from the statistical inference view 
[23] of the maximum-entropy principle where p*(p^) is seen as the 
least informative state compatible with the two-party marginals p^ 2 \ 
Since p and p*(p ^) are equal on every two-party subsystem of ABC, 
any discrepancy between them reveals additional information shared 
by p which can not be observed on any two-party subsystem. This 
information is called irreducible three-party correlation in [30] and is 
quantified by C 8 (p). 

We remark that information is seen as a constraint on our beliefs 
in the context of inference [10] as opposed to the language usage in 
coding theory where information is a measure of unpredictability which 
is quantified, in the case of quantum information sources [8], in terms 
of von Neumann entropy. 

The main point regarding continuity of C 8 is provided in Sec. 5.6 in 
[45]. Lemma 5.15(2) and Lemma 4.5 in [45] prove for all a G L(u) that 

p* : L( u) — > M. 8 is continuous at a 

(6.4) $ 

U 3 : Ai 8 —>• M is continuous on E|_ A/ ) 8 (a). 

The equivalence (6.4) allows us to apply convex geometric methods 
to detect discontinuities of C 3 . In what follows, inspired by a model 
presented in Example 6 in [13], we prove existence of a discontinuity. 
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Example 6.1. We consider the two-local Hamiltonians 

H 0 := 1 2 <g> <r 3 g) cr 3 + a 3 <g) 1 2 g) <r 3 + a 3 g) <r 3 g) 1 2 
Hi := on g 1 2 g) 1 2 + 1 2 g oy g 1 2 + 1 2 g 1 2 g <Ji 


where 0 \ and <r 3 are Pauli matrices and for e > 0 we take 


H{e) := H 0 + eH v 

The maximal eigenvalue A(e) := 1 + e + 2\/l — e + e 2 of H(e) is non¬ 
degenerate. The positive number s = s(e) := (A(e) — 3)/3e which 
goes to zero for e —* 0 allows us to write w(e) := (l,s,s,s,s,s,s,l) 
for the corresponding eigenvector which defines the pure state p(e) : = 
w(e)w(e)*/(2 + 6s 2 ) with AB-marginal 


p(£W=|(|00>(00|+ 111)011) + ^ 


—2s 2 s 2 +s s 2 +s 
* 2 s 2 2s 2 


2 s 2 


2s \ 
s 2 +s 1 
s 2 +s ) 
-2s 2 / 


G Ad 4 . 


The symmetry of H(e) shows p(e) (2) = (p(e) A B, p(c)ab, p(c)ab)- Using 
the identification (6.2) we note from Fact 5.1(2) that p(e)*U is an ex¬ 
posed point of the convex support L( u) where u = (ui,. . . ,u r ) spans 
the space of two-local Hamiltonians. The limit p(O)U) of p(e)^ for 
e —> 0 is the mid-point of the segment between the distinct marginals 
|000)(000| (2 ) and |111)(111|U) and therefore p(O)U) is not an extremal 
point of L(u). 

Out of curiosity we mention that p(e) A B has rank two for 0 < s < 1 
(non-zero eigenvalues an d ^jiyyrjy)- Although the extremal 

points of Ad 4 are rank-one states this does not contradict the fact that 
p(f) <2) is an exposed point of L(u) because L(u) C Ad 4 x Ad 4 x Ad 3 . 

Let us turn to the irreducible correlation and its continuity. We 
have just seen that p( 0 ) (2 ' ) is a limit of exposed points but not an 
extremal point itself. Hence Remark 5.7(1) shows that the maximum- 
entropy inference p* is discontinuous at p( 0 )*A and (6.4) proves that 
the irreducible correlation U 3 is discontinuous at some point in the fiber 
{p G Ad s | pW = p(O)U)}. 

This abstract existence result of a discontinuity is confirmed by the 
discontinuity of U 3 at the GHZ state ^(|000) + |111)) which follows 
directly from [30] and which we explain in detail. The projection 


p ■= | 000 )( 000 | + | 111 )( 111 | 


is the spectral projection corresponding to the maximal eigenvalue of 
H 0 . So Fact 5.1(2) applied to u(A) = H 0 proves that Ad(pM 8 p) is the 
inverse image of an exposed face of L(u). Elements of the Bloch ball 
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M (pMgp) can be written in the form 

a(x, y,z):=^(p + x- (|000)(111| + |111)(000|) 

+ y ■ (-*|000)(111| + i|lll)(000|) 

+ z- (|000)(000| - |111)(111|)) 

where (x, y, z) G M 3 lies in the three-ball, that is x 2 + y 2 + z 2 < 1. The 
two-party marginal is 

a(x,y,z)W = |(1 + ^)(|000)(000| (2) + ±(1 - z)|lll)(lll| (2) . 

The maximum-entropy state 

P*(p(x,y, -) (2) ) = K 1 + ^)( I ooo) (ooo | + |(i - z)|iii)(ni| 

has von Neumann entropy S(p*(a(x, y, z)^)) = R(|( 1 + z)) where we 
use the function H(rj) := — rj\og(rj) — (1 — rj) log(l — rj), q G [0,1]. In 
particular, for pure states^ := a|000)+ / 5|lll), a,/3 G C, |or| 2 H-1/31- = 1 
we get 

c 3 {ii>r) = H( m 2 ) 

which is strictly positive unless |a| = 0 or |a| = 1. The irreducible 
three-party correlation C^iftift*) has the maximal value log(2) for states 
ipip* where ift = ^(|000) + e*^|lll)), (ft G [0, 2 tt). 

On the other hand, ift is approximated, for example, by vectors p : = 
a|000) + cos(7)/3|111) + sin(7)/?|001) for real 7 —> 0. For small 7 > 0 
the vector p is not a local unitary transform of a vector ift because 
the two-party marginals of pp* are not identical. This implies [30], as 
we have recalled in (6.1), that pp* is uniquely determined by its two- 
party marginals. Hence (pip* belongs to the maximum-entropy states 
p*{L( u)) and has zero irreducible three-party correlation. This proves 
discontinuity of C 3 at iftift* for all 0 < |a| < 1. Analogous discontinuity 
statements hold for unitary transforms of vectors ift. 

7. Multiply Generated Round Boundary Points 

In the sequel we study pairs of observables u 1 , u 2 £ M£, dGN, with 
E u : M d h — y M 2 for u = (ui,u 2 ) and where E u (A4d) = L(u) = W(u) is 
the numerical range (2.4). In this section we prove a sufficient condition 
for the openness of E|x d in terms of unique pre-images and simplicial 
points. We will see that this condition works well for matrix sizes 
d = 3,4, 5 but has a limited meaning for d > 6. 

In our earlier work [28] we call z G W{u\ + iu 2 ) C C singularly 
generated if x*(ui+iu 2 )x = z holds for exactly one linearly independent 
unit vector x G C d . Otherwise z is multiply generated. Since E|_A/t d is 



L. Rodman, I. M. Spitkovsky, A. Szkola, S. Weis 


19 


open at those points p G A id where E(p) is an interior point of the 
numerical range, a classification of boundary points is useful. A corner 
point is a point of W{u\,u 2 ) C M 2 which belongs to more than one 
supporting line of W{u\,u 2 ). A flat boundary portion is a non-trivial 
line segment lying in the boundary of W(ui,u 2 ). A boundary point 
of W{u\,u 2 ) which is no corner point and which does not belong to 
the relative interior of a flat boundary portion is called round boundary 
point. 

Notice that every flat boundary portion is a sub-segment of a one¬ 
dimensional face of W(u\,u 2 ). Round boundary points exist only if the 
numerical range has dimension two. 

Lemma 7.1. If the dimension of W{u\,u 2 ) is two then every corner 
point is the intersection of two facets of W ( 111 , 1 / 2 ) ■ Without dimension 
restrictions, every point of W[u\, i/ 2 ) is either a round boundary point 
or a simplicial point but not both. 

Proof: Every corner point a G W{u\ + iu 2 ) is a normal splitting 
eigenvalue , see Sec. 13 in [24], that is u\ + iu 2 is unitarily equivalent 
to a block diagonal matrix 


a 



B 


with zeros on the off-diagonal. Since W{u\ + iu 2 ) is the convex hull of 
W ( B) and of a the first statement follows by induction. Going through 
the above classification of boundary points, the second assertion now 
follows easily. □ 

Remark 7.2 (Lattice theoretical proof of Lemma 7.1). The numerical 
range W(ui,u 2 ) is the convex dual of an affine section of Aid [21, 
42, 20]. Since all faces of this affine section are exposed, the lattice 
isomorphism (2) in [43] shows that all non-empty faces of normal cones 
of W(ui,u 2 ) are normal cones. In particular, if a is a corner point 
then the two boundary rays of its normal cone are normal cones of 
W( Ul ,u 2 ). The lattice isomorphism (1) in [43] now shows that a is 
an extremal point of two distinct one-dimensional faces of W{u\,u 2 ) 
which proves the claim. 

Theorem 7.3. The map E u |^ d is open on Aid except possibly at those 
states p G Add where E u (p) is a multiply generated round boundary 
point. 

Proof: If a G W(ui,u 2 ) is singularly generated then by Thm. 5 in 
[12] the fiber E|]^ (a) is a singleton. It is easy to prove and well-known 
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in the theory of multi-valued maps [15] that E| M d is open at singleton 
fibers. If a is a simplicial point then Prop. 3.2 proves that E| M d is open 
on E|^ (a). Otherwise, if a is not a simplicial point, then the second 
statement of Lemma 7.1 shows that a is a round boundary point. □ 

The converse of Thm. 7.3 does not hold for d > 4, that is, E|_Aa d may 
be open on fibers of multiply generated round boundary points. 

Example 7.4. An example inspired by Thm. 4.4 in [28] is 

. f 0 2 1 ro 2 1 „ [0 2 ' 

(7.1) U\ + iu 2 00 ® 00 I 2 ® 00 ’ 

where the numerical range is the unit disk and all points on the unit 
circle S 1 are multiply generated round boundary points. Nevertheless 
E|_yn 4 : Ad 4 —> M 2 is open. A short computation with the unitary 

vo := I 2 ® [cos(|)l 2 - i sin(|)<T 3 ] G M 2 , 6 G M, 

and the Pauli matrix cr 3 shows 

w cos(d) — sin(d) . . . 

nvopve) = [ sin( ^ cos( ^J-EW. p€M t . 

The state space Ad 4 is partitioned into the orbits of the group of uni- 
taries {vg \ 6 e [0,27r)}. The orbit O(p) through any state p with 
E(p) G S 1 is homeomorphic to S 1 under E. Hence E|o( p ) is open at p 
and a fortiori E restricted to E|) V | 4 (S' 1 ) is open at p. Now Thm. 4.3 
proves that E|_a ^ 4 is open at p. The openness of E |» 4 at all other states 
of M 4 follows from Prop. 3.2 applied to the interior of the unit disk. 

To capture the non-generic direct sum form (7.1) we introduce a 
definition from [26]. A matrix A G M d is unitarily reducible if A is 
unitarily equivalent to a matrix in block-diagonal form with at least 
two proper blocks. Otherwise A is unitarily irreducible. 

Unitarily irreducible matrices have at most d — 3 multiply-generated 
round boundary points when d — 3,4,5. See Thm. 3.2 in [28] for d = 3 
and notice that every round boundary point is an extremal point. For 
d — 4,5 see Thms. 4.1 and 5.7 in [28]. Since E|_a 4 2 is open by Fact 2.6 
and since Adi is a singleton, Thm. 7.3 implies the following. 

Corollary 7.5. Let d < 5 and let u 1 + m 2 be unitarily irreducible. 
Then there are at most max{0, d — 3} points a o/1U(ui,u 2 ) such that 
E|x d is not open on the fiber E| A | d (a). 

The maximum-entropy inference p* : IU(-Ui,w 2 ) —y A4 d is indeed 
discontinuous at the multiply generated round boundary point(s) of 
W(ui,u 2 ), if there are any, when d — 4 or 5 and when U\ + iu 2 is 
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unitarily irreducible. This happens because these points are isolated. 
So, approximating a multiply generated round boundary point a by 
singularly generated extremal points cu we observe that p*(cti ) has rank 
one while the rank of p*(a) is at least two. This proves discontinuity 
of p* at a. 

ft is known for matrix size d > 6 and irreducible matrix U\ + iu 2 that 
W{u\,u 2 ) may have infinitely many multiply generated round bound¬ 
ary points. As was noticed earlier in [28], Sec. 6, this is closely con¬ 
nected to the failure of Kippenhahn’s conjecture. On the other hand, 
one of us has recently shown [47] that for all d G N the maximum- 
entropy inference p* : W{u\,u 2 ) —* Aid has at most finitely many 
discontinuities. 


8. Numerical Range of 3 x 3 Matrices 

In the numerical range approach, based on [24, 26, 36], we character¬ 
ize points of discontinuity of p*(a) depending on the type of shape of 
the numerical range of 3 x 3 matrices: unitarily irreducible — ovular, 
ellipse, with a flat portion on the boundary, and unitarily reducible - 
triangle, line segment, ellipse, and the convex hull of an ellipse and a 
point outside the ellipse. 

See Example 4.18 in [45] for the problem of openness of E|>f(.A) for 
unitarily reducible u\ + iu 2 and the C*-algebra A C M 3 generated by 
U\ + iu 2 . Here we consider the full algebra M :i and arbitrary 3x3 
matrices. 


Theorem 8.1. Let ui,u 2 € Mg. If E\ Ms is not open on the fiber 
E|^ 3 (A) of a point z G W(u\ + iu 2 ) then after reparametrization (2.5) 
the matrix u\ + iu 2 has the form 



and W (A) is the unit disk in C. The map E|^ 3 is open on Aiz with the 
exception of the fiber E|[^ (1) which is a three-dimensional ball where 
E|x 3 is only open at the pure state vivl for t'i := 4^(1,1, 0). 

Proof. If E|_a 4 3 is not open on E\fJ l3 (z) then Thm. 7.3 shows that z 
is a multiply generated round boundary point. Thm. 3.2 in [28] shows 
then that U\ + iu 2 is unitarily equivalent to a matrix B © [z] where B 
is a unitarily irreducible 2x2 matrix and z is a boundary point of the 
ellipse W{u\ + iu 2 ). This and (2.5) allow us to transform the matrix 
U\ + iu 2 into the above matrix A. Now W[u\ + iu 2 ) is the unit disk 
and 1 is a multiply generated round boundary point. By Thm. 3.2 in 
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[28] all points 1 on the unit circle S 1 are singularly generated, so 
the fiber E|/j } 3 (a) is a singleton by Thm. 5 in [12]. At singleton fibers 
E|x 3 is open. The map E|_a 4 3 is open on the hbers of all interior points 
of the unit disk, see Prop. 3.2, so it remains to examine the exceptional 
point 1 G S' 1 . 

Since for any a G S 1 \ {1} the fiber of ct is a singleton we have 
H(a) = v a v* where v a := M(l, a, 0). Choosing any po G E|^ 3 (l) and 
maximizing the quadratic form g Po (p ) := — (p — po,p — po), p G Ai 3 , 
gives H( 1) = po- The restriction H\$i is continuous at 1 if and only if 
Po = V\v\. This, by Thm. 4.3, is also the condition that H is continu¬ 
ous at 1 and, by Fact 2.2, the condition that E |^ 3 is open at po- □ 


The maximum-entropy states with respect to the transformed ob¬ 
servables Ui — G\ © 1 and u 2 = cr 2 © 0 in Thm. 8.1 are 


P*( 1) 


1/4 1/4 
1/4 1/4 


©[ 1 / 2 ] 


and 


P*(«) 


1/2 cm 1/2 
a • 1/2 1/2 


© [ 0 ], 


a | = 1 , a / 1 . 


This proves discontinuity of p* at a — 1. 

We conclude with a remark about several observables ui,...,u r G 
M 3 , r G N, and u = (u\,... ,u r ). 


Proposition 8.2. Let F be a face of L( u) = E(Ad 3 ) and assume 
0 < dim( F) < dim(L(u)). Then F is a Euclidean ball of dimension 
one, two or three, and all relative boundary points w G rb(F) have 
singleton fibers E|]/j 3 (w). 

Proof: Let w be an extremal point of F. Then {w} C F C L are 
faces of L , properly included into each other. Since inverse images of 
faces are faces, Fact 5.1(1) shows 

E|x! 3 (F) = M (P 1 M 3 PO, EI^VW) = M(p 2 M 3 p 2 ) 

for projections p\ © p 2 . By the assumption of strict dimension differ¬ 
ences we get 

3 = rank(l 3 ) > rank(p!) > rank(p 2 ) > 1. 

Thus rank(pi) = 2 holds and so E|j^ (F) is a copy of the three- 
dimensional Euclidean Bloch ball. Now the claim follows easily. □ 


Prop. 8.2 implies the following (for several observables, r G Ft). 
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Corollary 8.3. Let F be a facet of L( u) = E(AL 3 ). Then E|^ 3 is 
open on E|^ 3 (F). 

An example where Coro. 8.3 gives new insights beyond Coro. 4.4 is 
the convex hull of the Steiner Roman surface, depicted in Fig. 10 in 
[4], which is a linear image of AI 3 and which has four disk facets. 


9. Appendix 

We recall a reduction of the state space in terms of the algebra of 
observables. The proof is from Sec. 3.4 in [42] and is reproduced here 
in a simpler setting. 

For d G N let us consider the parametrization of the hermitian pencil 

u(9) := 9\Ui + • • • + 9 r u r , 9 G IF 

for observables u 1 ,.. ., u r G M d . Let M(u) denote the real or complex 
*-algebra generated by the d x d identity matrix L d and by u±,... ,u r . 
Recall that Minkowski’s theorem asserts that a convex body is the 
convex hull of its extremal points. The statement of Straszewicz’s 
theorem is that the closure of exposed points of a convex body contains 
all its extremal points. See e.g. Coro. 1.4.5 and Thm. 1.4.7 in [37] for 
these statements. 

Lemma 9.1. Any spectral projection of any matrix in the hermitian 
pencil (u(0) | 9 G IF} belongs to M(u). The convex support of u is 
L(u) = E(M(M(u))). 

Proof: Any spectral projection of u(0) for 9 G IF is a real polynomial 
in one variable evaluated at u (9). Therefore the spectral projection 
belongs to the algebra M(u). 

Let a be an exposed point of L(u). Then Fact 5.1(2) shows that there 
is 0 G IF such that the spectral projection p of u(0), corresponding to 
the largest eigenvalue of of u(9), yields 

ei : m >) = MpM dP ). 

Since p/ tr(p) belongs to the algebra A := M(u), the image of the 
state space A4(A) under E covers all exposed points of L(u). Now 
Straszewicz’s theorem implies that A 4(A) covers all extremal points of 
L(u), and thus by Minkowski’s theorem the whole L(u). The converse 
inclusion is trivial because A C M d holds. □ 
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